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ABSTRACT 

We derive a relation for the steepening of blazar 7-ray spectra between the 
multi-GeV Fermi energy range and the TeV energy range observed by atmo- 
spheric Cerenkov telescopes. The change in spectral index is produced by two ef- 
fects: (1) an intrinsic steepening, independent of redshift, owing to the properties 
of emission and absorption in the source, and (2) a redshift-dependent steepening 
produced by intergalactic pair production interactions of blazar 7-rays with low 
energy photons of the "intergalactic background light" (IBL). Given this rela- 
tion, with good enough data on the mean 7-ray SED of TeV Selected BL Lacs, 
the redshift evolution of the IBL can, in principle, be determined independently 
of stellar evolution models. We apply our relation to the results of new Fermi 
observations of TeV selected blazars. 

Subject headings: Gamma-rays: general - blazars 

1. Introduction 

Stecker & Scully (2006) (SS06) derived a simple analytic expression for the change 
in spectral index of a TeV 7-ray source in the redshift range between 0.05 and 0.4. They 
showed that the change in the spectral index caused by intergalactic absorption is given by 
an approximately linear relation in redshift, i.e., Ar a ~ C + Dz. 

The purpose of this letter is to generalze this relation by including the effect of intrinsic 
steepening in the source spectra between the mutli-GeV energy range observed by Fermi 
and the TeV range observed by atmospheric Cherenkov telescopes. Our general result is 
roughly independent of the specific model of the intergalactic background light (IBL) used, 
because it only depends on the shape of the average galaxy spectral energy distribution 
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(SED) on the near IR side of the starlight peak that determines the absorption in the TeV 
energy range. 

We compare our relation for the specific baseline and fast evolution models of Stecker, 
Malkan & Scully (2006) (SMS06) with the results of recent Fermi observations of 13 TeV 
selected BL Lac AGN. We also show how it can be used to independently determine the 
redshift evolution of the IBL. 

2. Steepening by Absorption 

In order to determine the effect of intergalactic absorption, we use the results of SS06 
demonstrating that t(E 1 , z) can be fitted to an approximately logarithmic function in E 1 
in the energy range 0.2 TeV < E 7 < 2 TeV and one which is linear on z over the range 
0.05 < z < 0.4. It is important to note that our linear fit to the z dependence is both 
qualitatively and quantitatively different from the linear dependence on redshift which would 
be obtained for small redshifts z « 1 and which simply comes from the fact that for 
small z where luminosity evolution is unimportant and where r oc d, with the distance 
d ~ cz/H oc z. Our quantitative fit for the higher redshift range 0.05 < z < 0.4 comes 
from the more complex calculations based on the models of SMS. For this reason, the linear 
fits in SS06 are not simply proportional to redshift. 

SS06 found that t(E 7 , z) is well approximated by 

t(E 7 , z) = (A + Bz) + (C + Dz) ln[£ 7 /(l TeV)], (1) 

where A, B, C and D are constants. This expression holds over the energy and redshift 
ranges given above. The energy range of validity is the energy range to which the 
atmospheric Cerenkov telescopes are sensitive. 
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Following SS06, we assume an intrinsic source spectrum that can be approximated by 
a power law of the form $ SjTe y(-E 7 ) ~ kE 1 a ' TeV over a limited energy range. Then the 
spectrum that will be observed at the Earth following intergalactic absorption will be of the 
power-law form 



This can be compared with the empirically observed TeV spectra that are usually 
presented in the literature to be good fits to power-laws. The observed spectral index, 
F ,TeV, will then be given by 



where the intrinsic spectral index of the source is steepened by intergalactic absorption by 
an amount AT a (z) = C + Dz. 

On the other hand, in the multi-GeV range over which the Fermi Large Area Telescope 
(LAT) is sensitive, we expect essentially no steepening from absorption over the redshift 
range of validity, i.e., AT ajGeV (z) ~ 0. 

The parameters C, and D obtained by fitting the optical depths derived for the fast 
evolution (FE) and baseline (B) models of SMS06 are given in Table 1. 



<5>o,Tev(E,) = ke^ A+B ^E-^^ v+c+Dz \ 



(2) 



r ,TeV — T SyT eV + ^a,Tev{z) 



(3) 



3. Intrinsic Steepening 



The importance of synchrotron self Compton emission in producing high energy 7-rays 
in astrophysical sources was pointed out by Rees (1967). Particular applications to the 
Crab and other sources were discussed by Gould (1965) and Rieke & Weekes (1969). Today, 
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given present observational studies of TeV-selected BL Lac AGN, it is generally accepted 
that the synchrotron self Compton mechanism is the primary emission mechanism for 
producing 7-rays in the TeV energy range in TeV-selected BL Lac AGN. 

The 7-ray spectrum from Compton interactions is a smoothy varying one (Blumenthal 
& Gould 1970). Good empirical fits the the SEDs of the low energy peaks in blazars were 
obtained by using parabolic functions on log-log plots (Landau, et al. 1986, Sambruna, 
Maraschi & Urry 1996). In fact, these considerations led to the prediction of candidate 
BL Lac objects to be found by atmospheric Cerenkov telescopes in the TeV energy range 
(Stecker, de Jager & Salamon 1996; Costamante & Ghisellini 2002). 

We assume such an approximation here, i.e.. 



e2 ^e = /(a;) ' (4) 

where x = \ogE with the simplification that E is in dimensionless fiducial units (i.e., E/Ef), 
and 



/(*) = z %r L + b. (5) 

Here, A = log < W > is a constant parameter based on the ensemble averaged 
width of the parabola in log E space based on a set of SEDs of TeV selected BL Lacs and 
x = log Eq where E is the energy at which the Compton peak of the average BL Lac SED 
is a maximum. We note that for some BL Lacs electron scattering in the Klein-Nishina 
regime will distort our parabolic assumption on the high energy end of their SEDs, possibly 
reducing log < W > by a small amount (Tavecchio, Maraschi & Ghisellini 1998). However, 
this does not significantly affect our formalism. 

We find that the change in the slope of the intrinsic source spectrum occurring between 
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some average observed energy in the GeV range < -EceV > and some average observed 
energy in the TeV range < E Tc y > is approximately given by a constant, 

log < E TeV > - log < E GeV > 

A = K (6) 

Thus, the ensemble averaged intrinsic BL Lac source spectrum in going from the GeV 
range to the TeV range is independent of redshift within the errors in the observational 
data. We will therefore approximate the average intrinsic steepening over an observed set 
of TeV selected blazars by 

F s ,TeV — ^s,GeV = AT s = K. (7) 

4. The Relation for the Total Steepening between GeV and TeV Energies 

The total steepening expected betwen the GeV and TeV energy ranges is then given 
by the relation 

Ar = Ar s + Ar a = (c + k) + d z (8) 

The parameters C, and D obtained by fitting the optical depths derived for the fast 
evolution (FE) and baseline (B) models of SMS06 are given in Table 1. The parameter K is 
then derived by performing a x 2 fit to the observed steepening data obtained by the Fermi 
collaboration (Abdo, et al. 2009) in order to find (C + K). The parameter K is also given 
in Table 1 for the FE and B models. The data on the TeV selected BL Lacs from Abdo, 
et al. (2009) are shown in figures 1 and 2, along with the best-fit linear relations for the 
models indicated. It can be seen that, given the present limited data set and large error 
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bars, one cannot uniquely determine the parameters C, D and K. However, in principle, 
with a good enough data set, one could determine K uniquely from equation Then, 
since K ^> C, one could use equation (jSj) to determine the parameter D. This will then 
give a determination of the redshift evolution of the IBL independently of models of the 
star formation rate. 

Figures [I] and [2] show the fits of the parameters given in Table 1 to the linear 
dependence in and redshift given by equation flS}. 

5. Conclusions 

We have derived a simple analytic approximation for determining the steepening in the 
spectra of TeV selected BL Lac AGN and compared our results with recent observational 
data on 13 TeV selected BL Lac AGN. Our relation is in excellent agreement with the 
observational data and provides a framework for understanding and interpreting both these 
and observational data. 

SS06 have shown that the effect of intergalactic absorption on the spectra of AGN 
in the energy range 0.2 TeV < E 1 < 2 TeV and the redshift range 0.05 < z < 0.4 is a 
simple power-law to power-law steepening. Absorption in this energy range is primarily 



Table 1. Steepening Parameters 



Evolution Model 


C 


D 


K 


Fast Evolution 


-0.0972 


10.6 


0.427 


Baseline 


-0.0675 


7.99 


0.716 
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Fig. 1. — The fits obtained for the linear functions C + Dz (dashed line) and C + Dz + K 
(solid line) shown for the SMS06 baseline model as descibed in the text. These are fit to the 
data on 13 BL Lacs given by Abdo et al. (2009). 

produced by interactions with near infrared photons from on the low energy sides of the 
starlight peaks in galaxy SEDs. The shape of the resulting peak in the intergalactic SED 
produces an approximately logarithmic energy dependence for the function r(E y ). This 
energy dependence should be roughly the same for all models of the IBL. Therefore, our 
general result of a linear Ar = (C + K) + Dz relationship is roughly independent of the 
specific IBL model used because it only depends on the shape of the average galaxy SED 
on the near IR side of the IBL starlight peak. The parameters C, D and K will be different 
for the different models since the absolute value of r varies from one IBL model to another. 

Given our derived relation between Ar and redshift, and with good enough data on 
the mean 7-ray SED of TeV Selected BL Lacs used to determine K, the redshift evolution 
of the IBL can, in principle, be determined independently of stellar evolution models. 




1 0.00 0.05 0.10 0.15 0.20 0.25 0.30 

Redshift z 

Fig. 2. — The fits obtained for the linear functions C + Dz (dashed line) and C + Dz + K 
(solid line) shown for the SMS06 fast evolution model as descibed in the text. These are fit 
to the data on 13 BL Lacs given by Abdo et al. (2009). 
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